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Abstract. A Spin filtering device through quantum spin interference is addressed, 
in two dimensions, in a GaAs/AlGaAs electron gas that has both Rashba and 
Dresselhaus spin-orbit couplings and an applied external magnetic field. We propose 
an experimentally feasible electronic Mach Zehnder Interferometer and derive a map, 
in parameter space, that determines perfect spin filtering conditions. We find two 
broad spin filtering regimes, one where filtering is achieved in the original incoming 
quantization basis, that takes advantage of the purely non-Abelian nature of spin 
rotations, and the other, where one needs a tilted preferential axis to observe 
the polarized output spinor. Both solutions apply for arbitrary incoming electron 
polarization and energy, and are only limited in output amplitude by the randomness 
of the incoming spinor state. A full account of beam splitter and mirror effects on spin 
renders solutions only on the tilted basis, but encompasses a broad range of filtering 
conditions. 
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1. Introduction 

The Rashba and Dresselhaus SO interactions arise in materials which lack either 
structural or bulk inversion symmetry, respectively [1] [21 |3]. These two kinds of 
interactions have recently been given a great deal of attention due to their potential 
role in the generation and manipulation of spin polarized currents, spin filters|ll[5l|6l[7], 
spin accumulation [8], and spin optics [9]. 

A reformulation of the spin-orbit coupling Hamiltonian in terms of non-Abelian 
gauge fields [To] was explicitly given in ref. [11] |12l [131 HI] where the SO interaction 
is presented as a SU{2) x U{1) gauge theory. As the Yang-Mills gauge theory is well 
understood and is the underpinning of well established theory, enormous insight can 
be brought upon new problems. Such gauge point of view, in more general terms, 
has been known for some time|15l [T6l [T7] . This formulation is very revealing, since 
the consistent gauge structure of the theory becomes obvious and the physics of spin 
currents, persistent currents and color diamagnetism [18] can be understood in a manner 
analogous to the well known U{1) gauge theories. A consistent SU{2) x U{1) gauge 
approach was presented in reference [131 [E] where it was found that for the Pauli 
type Hamiltonians (including Rashba and 2 dimensional reductions of the Dresselhaus 
Hamiltonian), Gauge Symmetry Breaking (GSB) is necessarily built into the theory and 
leads to vanishing of the spin conductivity in constant electric fields [ll]. In addition, 
the Yang Mills interpretation of the Rashba and Dresselhaus SO interactions renders the 
associated gauge fields real, with topological consequences analogous to the Aharonov 
Casher effect [131 Hi] . 

Recent proposals were recently reported for the construction of perfect spin filters 
based on active Rashba spin orbit media [6], ballistic spin interferometers [T9] and the 
analysis of the persistent spin helix [20 [ [7], where the Yang Mills gauge point of view 
is advantageous. Here we readdress the problem of spin filtering by interferometry in 
a quasi two dimensional system, and make connection to an experimentally feasible 
test of these ideas through an electronic Mach Zehnder interferometer (MZI) within 
Rashba and Dresselhaus media. Recent proposals contemplating this setup as an 
spin intereference device include quantum logic gates[2T], bit controlled Stern-Gerlach 
devices[5] and tunable entanglement [22] . Our analysis, within this setup, enables us to 
obtain exact conditions for spin filtering which can be achieved by tuning appropriate 
experimental parameters. Such conditions for spin filtering greatly generalize previous 
special situations where the spin polarization is a conserved quantity [23], and show new 
possibilities for spin filtering beyond previous approximate treatments. 

The structure of the paper is as follows. First we consider the Hamiltonian with 
both Rashba and Dresselhaus interactions for a two dimensional electron gas (2DEG) 
including a magnetic flux described by a [/(I) gauge field. Following the approach given 
in ref. [7j, we rewrite the Rashba and Dresselhaus contributions in terms of a Yang 
Mills gauge field and review how this approach leads to the introduction of a GSB term 
analogous to the Proca term for the Maxwell field. Then, we propose an interference 



Perfect spin filtering device 



3 



setup in the form of an electronic MZI where the electron's spin transport is modulated 
due to the presence of Rashba and Dresselhaus active media. We derive the conditions 
for perfect spin filtering that are applicable independently of the incoming spin state 
and the full energy range of the injected electrons. Finally, we give some concluding 
remarks. 



2. Spin-Orbit scattering for two dimensional electron gas 

We consider a two dimensional system consisting of non interacting electrons subject 
to both Rashba and Dresselhaus spin orbit interactions. In addition, one can apply an 
external transverse magnetic flux $5 described by a U{1) gauge vector potential A. 
Two recent works have shown how to measure and control the Rashba and Dresselhaus 
parameters using gate voltages in two dimensional GaAs/AlGaAs electron gas[2ll[25]. It 
is striking that one can achieve SO magnetic fields of 2-3 mT. The SO physics beautifully 
follows an extended weak localization theory that allows for a detailed access to the 
material parameters. 

One can address the two dimensional GaAs/AlGaAs electron gas by a single particle 
Hamiltonian including the previously described couplings by 

TT2 fc^ 

^ = ^ + ^- "(n-^' - n,a^) - /3(n,a^ - U.a^ + -^o^ (1) 

where 11 = p + eA, — e and m* are the electron's charge and effective mass, V 
a substrate lattice potential that can be assumed periodic, cr is a vector of Pauli 
matrices, and a and /3 are material-dependent parameters characterizing the Rashba 
and Dresselhaus interactions, respectively. The last term is the Zeeman energy. The 
term linear in k describing the Dresselhaus interaction results from averaging a cubic 
in k contribution (for the bulk) in the confining direction and neglecting other cubic 
terms in the strong lateral confinement situation |26]. In the rest of this work we ignore 
the effect of the Zeeman term in the limit of small magnetic fields (a few flux quanta 
through a 200 x 200/im^ area) such that the spin orbit energy is much larger than the 
Zeeman energy [27j. According to measured parameters in ref. [21] the SO energy for 
an GaAs/AlGaAs electron gas is 5 orders of magnitude greater than the Zeeman energy 
for the proposed field strengths. This way the external magnetic field results in strong 
phase effects through the vector potential but no appreciable precession occurs due to 
the Zeeman term. Nevertheless, we will see that there are spin filtering scenarios for the 
device even for zero external magnetic field. 

Following [6], [7] , we introduce a spin dependent (non-Abelian) gauge field W whose 
components are given by 

= (/3r^ - arJ')x + (ar^ - /3r^)y, (2) 

m* 

with t"" = cr"/2, and g/h is the SU{2) coupling constant. Using this gauge field we can 
rewrite equation ([1]), having ignored the Zeeman contribution, in the form 

(p + eA + gWr"? , g'^W ■ W 
H = h cAq . (3j 

2m* 8m* 
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Figure 1. Sketch of the electronic Mach Zenhder interferometer setup. The arms of 
the square are made of active SO Rashba and Dresselhaus media. The beam splitters 
are implemented through two Quantum Point Contacts (QPCs). There is a magnetic 
flux $B through the square. 



The first term describes the total kinetic energy taking into account the contribution 
from the regular vector potential due to an external magnetic field and the non Abelian 
gauge field. The second term is the background lattice potential whereas the third 
term represents a gauge symmetry breaking contribution similar to the field originally 
discussed in references [131 El l2H] responsible for rendering the spin currents physical. 

3. Electronic Mach Zehnder spin interferometer 

A device configuration that allows us to address the problem of spin filtering in a gauge 
independent |14] manner is the Mach Zehnder Interferometer (MZI). The setup for an 
MZI is sketched in (figure [1]). Here we are interested in determining the resulting 
amplitude ^j/u^ at detector Di, with i = 1,2 and to find the conditions for perfect spin 
filtering [6] at either detector. There is an interesting issue that must be discussed 
regarding spin 1/2 filtering. If the state at the input is a pure state spinor of spin 1/2, 
the electron is polarized on some indeterminate axis, in principle random, coming from 
the Fermi sea of the input conductor. If one could find this axis for every electron 
extracted then one would have a perfect spin filter for each electron. Nevertheless the 
resulting current is unpolarized. We thus define the spin filter as one acting on any 
entering (pure state) polarization and returning a polarized state along a definite axis. 
This approach will serve to build a polarized spin current. 

The relevant processes within the interferometer are described as follows (see figure 
[I]): Single electrons are assumed to be extracted from the Fermi sea as pure states 
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\l/o = (^-) • The electrons then pass through a beam sphtter that can be implemented by 
a combination of Quantum Point Contacts[29] the first of which we label QPCi described 
by a 4 X 4 scattering matrix Si that mixes spin orientations on perpendicular reflection, 
while it is diagonal for direct (no change in direction) transmission |30] . Mixing of spin 
orientations occurs at all reflections (including mirrors) due to changes in direction of 
the electron k vector within spin-orbit active media that changes the orientation of 
the implied wavevector-dependent magnetic field. Furthermore, as we consider both 
Rashba and Dresselhaus interactions, we need to derive general reflection conditions at 
the beam splitters and mirrors. In reference [30], this was done for Rashba assuming 
that small enough spin-orbit strength would yield only a small divergence of the reflected 
spin states in a k dependent basis. Surprisingly, when only the Rashba interaction is 
involved, the reflection matrix depends only on the incident angle and the reflection 
coefficient. On the other hand, if both Dresselhaus and Rashba are included, this is 
no longer true, and except for special angles of incidence, the reflection matrix depends 
on both Rashba and Dresselhaus strengths. The general reflection matrices are derived 
in the appendix. In this paper we will take the limit of 7r/4 reflections, that leads to 
simple, spin-orbit independent matrix elements. 

The resulting beams follow path I {II) that consists of a first horizontal Cj (vertical 
Cii) arm made of Rashba-Dresselhaus medium whose length is Lj (Ljj). The electrons 
are then specularly reflected from an ideal mirror Mi (M2), that also mixes spin 
directions, followed by a vertical £'/ (horizontal C'jj) arm of length Lj (Ljj) of the same 
material. The mirrors can be implemented as a simplified version of the beam splitters 
of reference [29j. Then the electrons pass through a second QPC (QPC2) described by 
the corresponding S-Matrix 5*2. Finally, two electron beams are collected at detector 
Di {i = 1,2), and we have = "ifji + "^n^i, where, \E'/^j (^//,i) is the corresponding 
transferred spinor through the it/i-arm. These amplitudes can be written in terms of 
the injected spinor \E'o as \E'd. = Uoi'^o, where the 2x2 matrices Ud, (generalized 
comparator operators [31]) are given by 

" i 



= (^2) exp 



d\-{p - eA - gWr") (r^) exp 



{r2i) exp 



D2 



[r2r) exp 



(^2) exp 



^ / d\-(p-eA- 

^ Jc'n 
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Such operators applied to the initial state do not change the energy expectation 
value. The transmission and reflection matrices regarding both Rashba and Dresselhaus 
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interactions, for it/ 4 incidence angle, are given by 




±irj 



) 



(5) 



where the subscripts j correspond to the beam sphtter index (see figure [T]) and r, I 
(corresponding to +, — in the non diagonal matrix elements, respectively) encode 
whether the electron current is refiected counter-clockwise (/) or clockwise (r). and 
tj are the refiection and transmission coefficients for the j— th beam splitter, while for 
the mirrors, the reffection coefficients are equal to 1. Note that Un- is not a unitary 
operator. The normalization condition |^I/diP + I^DjP — ^ ^^e total probability 
at the detectors require that Uj^JJc-^ + U\)jUD2 = \ the unit matrix. This simply 
means that the amplitudes received at the detectors do not interfere. The arms of the 
interferometer can be built from gate defined quasi one dimensional paths implemented 
on a 2DEG, where all transport is kept within one of the available transverse modes. 
The scattering length is assumed to be long enough, so that phase relations can be 
accurately described by the path lengths and the spin-orbit strengths as in the Datta 
Das[32] switch arrangement. 

4. Results: Spin diagonal mirrors and beam splitters 

In this section we consider a simplified version of the filtering device where beam splitters 
and mirrors are considered diagonal matrices or scalars. Although this approximation 
does not contemplate the matrix nature of the refiections we will obtain a simple scenario 
for the filtering properties of the device. The full problem will be treated below where 
essentially the same quahtative results are obtained. 

If the electric field E is uniform and static, the operators p — eA and (^iVV^t" 
commute. Thus, we can separate the orbital from the internal translation operators. 
For simplicity we will assume a square interferometer, thus Lj = Lu = L. Otherwise 
there are no restrictions or approximations related to the dimensions of the arms of 
the interferometer. As in Chen and Chang [7J we will make the discussion general by 
treating both the Rashba and Dresselhaus spin-orbit coupling on equal footing. 

Concerning the orbital contribution, it is easy to see that this will consist of a 
global phase exp[p ■ (Li + L2)] which we can drop, and a relative U{1) phase ipB which 
arises from the noncommutation of p and A. Using the definition for the magnetic 
fiux $B = BL^ and that for the fiux quantum 0o = h/e, the nontrivial orbital phase 
is written as 27iipB = 27t^b/4'o- On the other hand, the internal part gives rise to the 
SU{2) spin-dependent phase contribution. In order to simplify the resulting expressions, 
we introduce the adimensional variable 



that will be the crucial control parameter governing the SO interaction. Furthermore, we 
introduce the definitions 6 = tan~^(/3/Q;) along with the matrices 5"^ = cos 6*0"^ — sin6'cr^ 



(6) 
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and = sin^o"^ — cos6cry, such that (5"*)^ = 1, with 1 the identity matrix in spin 
space. After the previous considerations we can rewrite equation H] in the form 

Kd, = (t2)exp(-iAa^)(n)exp(-iAa2)(ti) + 

exp{27ri(pB){r2r) exp(-zAcr^)(rr) exp(-zAcr^)(ri/), 
= (r2r)exp(-iAa-^)(r/)exp(-iA5-^)(ti) + 

exp{2iTTipB)(t2) exp(— 'iAo"^)(rr) exp{—iAa^){rii). 

Due to the symmetry of these expressions (Udz is obtained from by the substitutions 
^"2 ^ ^2) we can focus on the first process, and obtain the second by making the necessary 
substitutions. Using the identity exp(±i7(T") = cos 7 11 ± ia'^ sin.'j, vahd also for our 
redefined a, the matrix W/)^ takes the form 

Uf)^ = titslcos^ A a - i sin A cos A (cr^ + a^) - ct^ct^ sin^ A] + 

rirse^^^^^ [cos^ Al-i sin A cos A(a^ + a^) - a^a^ sin^ A]. 

Now, we can easily determine that a^a'^ = sin 20 11 — ia^ cos 20 thus a'^a^ = sin 20 11 + 
icT^ cos 20 and (T^+a^ = (cos 0+sin 0)((t^— o"^). Substituting these results and rearranging 
the obtained expressions leads to 

Ujj^ = A+ [cos^ A - sin^ A sin 20] 11 + z sin A M, 

where we have introduced the traceless matrix M = A- sin A cos 20a^ — ,4+ cos A (cos + 
sin0)((T^ — 0"^) and A± = tit2 ± rir2e^*'^'^^ . The traceless condition simplifies the 
diagonalization of M, and the eigenvalues for are easily found to be 

A£i = A+[cos'^ A - sin^ A sin 20] ^ z sin A^^^ ^^^2 ^ ^^^2 20 + 2AI cos^ A(l + sin 20) (7) 

If we now define B± = tir2 ± rit2e^*'^'^^, the eigenvalues of the matrix are obtained 
from the previous result by making the substitution A± — )■ B± 

[cos^ A — sin^ Asin20] =1= i sin A-y/ sin^ A cos^ 20 + 2i3^ cos^ A(l + sin 20). (8) 

In order to get more insight into the nature of the conditions for perfect spin filtering 
we will specialize the previous expression to symmetric beam splitters i.e. ri = r2 = r, 
and ti = t2 = t. Within this case, we have A± = ± r^e^*'^'^^. Since we are interested 
in filtering one spin component, say the up component, we now proceed to determine 
the vanishing conditions of the corresponding eigenvalue A+\ 

From expressions ([71 (HI), these vanishing conditions can be found by either having 
cos A = or cos A 7^ (see also equation [6]). Although the former condition is 
mathematically only a particular case of the general solution, we distinguish it because 
the corresponding Un-^ becomes diagonal with respect to the original quantization axis, 
so we can speak of filtering along a non-tilted axis. Such a solution is also the simplest 
from the detection point of view since it involves the choice of a single quantization axis 
for the whole setup. The second condition (cos A 7^ 0) corresponds to finding a new axis 
where the up spin is filtered and we call such axis the tilted quantization axis. Note that 
both these filtering conditions (non-tilted and tilted) are independent of the polarization 
axis and the energy of the incoming state. We will comment further on this below. 
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4.1. Non-tilted filtering 

Let us first analyze the non-tilted situation. In this case the filtering condition requiring 
A;^^ = for all incoming energies (see equation [7]) , leads to the relation 

i(t^ - r^e^^'^'PB) 
tan2e = — f . 

Two 50 — 50 beam splitters for which r = i/ a/2, t = 1/ a/2, will then lead to the relation 
sin vry^B sin 26^ = cosvri/)^ cos 26', equivalent to the simple expression cos{7!'ipB + 26*) = 0, 
satisfied by the condition 

7i^B + 2e = {2n + l)^, (9) 

where n is an integer. Figure [2] depicts the relation between the spin-orbit parameters 
and the magnetic flux, for n,l = 0, necessary for perfect filtering of the up component 
in the original quantization axis. The spin-orbit parameters are in a reasonable range, 
as depicted in the figure, since for a GaAs heterostructure ha ~ 3.9 x 10~^^eV m[32]. 
h(3 ~ 2.4 X 10~^^eV m and h^/m*L ~ 1.7 x lO^^^eV m, assuming the arm of the 
interferometer ~ Ifim and an effective mass of m* = 0.046mo. These parameters yield 
< 6 in units of h/ {m*L). Note that our definition of a, P differs by a factor h to 
the standard definition (see equation [1]) . In reference |24] it is shown that gate control 
can vary a and /3 parameters by a factor of 6 by applying gate voltages in the hundreds 
of mV. 

The solutions are on a helix, as can be shown from the previous relations where 

h 



a 



V(2/ + l)7r/2cos[7r/4(2n + 1 - 2^^)], 



m*L 

(3 = V(2/ + l)7r/2sin[7r/4(2n + 1 - 2ipB)]. (10) 
m*L 

The integer n was defined in equation [9] while the second integer / is defined by the 
condition cos A = 0. 

The previous conditions, depicted in figure [21 do not tell us about the intensity 
of the signal received in detector Di i.e. the efficiency of the filter given an incident 
intensity. For this, one has to look back at the eigenvalues. While A:^^ = the amplitude 
of the outgoing polarized spinor at detector Di is given by 

^""^ " (a^'^^o") " Ge^-^^ cos(7r(^B - 2^)V'o") ' ^^^^ 

whose modulus squared is cos'^{n lpb — '^0)\'4^o\'^- Figure [3] shows a polar plot for the 
amplitude of the filtered signal (radius vector) as a function of the parameter designating 
the field flux (fB and the a, fi combination given by equation [10] for n = 0, 1 and / = 0. 
The figure shows that while filtering occurs for all the fluxes (given the appropriate 
values of a, /3) the amplitude can be zero, or very small, for some flux values i.e. in this 
case, the detector D2 gets most of the total amplitude. On the other hand, for some 
values of the flux, filtering can be very strong since the probability for a polarized spin 
can approach unity. The behavior of the second detector D2, while the first detector 




Figure 2. Perfect filtering for the non-tilted axis (original incoming basis). The 
plot shows the relation between a, /3 in units of h/{m*L), and ips that yields perfect 
polarization of the spin from an unpolarized input. The figure corresponds the values 
n,l = according to equation [TUl 




Figure 3. Filtering probability for the non-tilted solution of detector Di for n = 0, 
I ~ solid (blue) curve and n=l, 1=0 dashed (red) line. The radius vector depicted 
shows the filtered probability for the output spinor for one whole period in the 
parameters a, /3 as given in the figure[2j The position of the dashed vector corresponds 
to ifB = 0.25. The grey points represent "spin flipping" or opposite filtering solutions 
for detector D2. 
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sees a filtered signal, can be obtained through the eigenvalues of that detector having 
substituted the condition A:^^ = 0, namely 

ypj ^ -^in^B sin(7rv^B " 26). (12) 

It is obvious that the second detector D2 does not filter concomitantly with the Di in 
general. Furthermore, one can only find conditions for the second component to be zero 
(opposite filtering to detector Di) since the first component has modulus one. This 
takes us to the non-tilting spin flipped or opposite filtering solution at detector D2, only 
occurring while detector Di is filtering with maximal efficiency i.e. maximal polar radii 
in figure [3l 

The filtering amplitude is proportional to the projection of the incoming spinor 
(which has arbitrary weights onto the chosen quantization axis) to the surviving 
component at the output (see equation [TTj) . This means that for each arbitrary incident 
spinor from the Fermi sea one gets a filtering probability that depends on this projection. 
The resulting polarized current will thus have a random noise associated with this effect 
besides the contribution from shot noise. 

It is important to note that this solution does not appear in Abelian approximation 
(only exact in the case = and in one dimension) to the translation operator, where 
the SU{2) gauge vector operator has the same algebra as the U{1) gauge vector. The 
previous approximation was implemented in reference [7] by neglecting the commutator 
between components of the SU{2) gauge vector within a finite difference scheme. In 
this sense, the non-tilted case is an intrinsically non- Abelian scenario for spin filtering. 



4-2. Tilted filtering axis 

The tilted axis filtering scenario was discussed, within the tight-binding model, by 
Hatano, Shirasaki and Nakamura[6] when the Rashba coupling is present. In their 
approach, the interferometer involves an incoming lead and one outgoing lead, in 
contrast to our Mach-Zehnder configuration. The non- Abelian treatment is exact within 
their model, and requires a tilted outgoing axis to realize perfect spin filtering. 

For the Mach-Zehnder configuration, addressed here, the tilted axis solution (i.e. 
cos A 7^ 0), requires A^^ = 0, which implies 

A+ [cos^ A - sin^ A sin 261] = i sin A ^J Al sin^ A cos^ 29 + 2Al cos^ A(l + sin 29) . 
Squaring both sides and after some algebra one finds 

Al = sm^Acos^29{Al-Al). (13) 
Using the definitions for A±, and taking the square root, we reduce equation [13] to 

^2 ^ ^2^2in^B ^ 2rte^^'^s sin^ A cos 29. 



Employing the 50 — 50 mirror condition, we get after substitution 
sin TUfB = sin^ A cos 29. 



(14) 
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Figure 4. a) Perfect filtering by interference for the tilted axis. The plot shows the 
relation between a, (3 in units of h/{7n*L), and svcvklpb in a contourplot, the darker 
regions indicate larger values for the magnetic flux needed to yield perfect filtering, 
from an unpolarized input. Highlighted circles depict the zero flux solutions that yield 
perfect filtering, b) Perfect filtering probability for the tilted axis. The plot shows the 
relation between a, (3 in units of h/{ni*L), and the filtered intensity in a contourplot. 
The lighter regions indicate larger values for the intensity of filtering for the relation 
between parameters depicted in figure H^. Note that the circles evident from figure 
correspond to zero output amplitude. 



This is the relation between the spin-orbit parameters and the magnetic flux that 
leads to perfect filtering in the tilted axis. The solution is depicted in a contourplot in 
figure where the value of sin tci^b is represented in shades of gray as a function of 
a and /3. Each contour corresponds to a constant magnetic flux value and runs over 
the perfect filtering values of a and (3. The circular contour, depicted in the figure, 
corresponds to a = solution to equation [H] that leads to {m*L/K)^Ja^ + /3'^ = pir, 
for p integer. The figure depicts the solution for p = 1, 2, i.e. circles in units of h/ {m*L). 

In order to see if the filter is actually working, we must address the filtered 
amplitudes by looking to the second eigenvalue at detector Di. For the filtering condition 

Af?i = -2ze'^'^^ sin t^^^b [cos^ A - sin^ A sin 26] . (15) 

Substituting equation [14] in this expression and computing the modulus squared of the 
eigenvalue, we determine the strength of the filtered output, as was done in equation [TTl 
We have depicted the analytical solution for a range of values of a, /3 in the contour plot 
of figure Hb. The darkest shade corresponds to zero amplitude, and as the shade lightens 
the probability is higher for the filtered output. We note that the filtering solutions for 
the circular contours in figure and the lines a = ±/3 have zero amplitude. Such zero 
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Figure 5. Detector D2 output while Di filters out the up spin component (spin 
down polarization). The plus (minus) zones represent the regions where only the up 
spin (down) survives at the D2 detector. Note that either one or the other is filtered. 
The white regions represent no output in the detector and correspond to the localized 
phase. On can have either up or down spin filtering in D2 while up spin is filtered out 
in Di. 

amplitude solutions correspond to those of "localized solutions" of Cheng and ChangjT] 
where there is no filtered output. Behavior of detector D2, while Di is filtering out 
the spin up component (spin down polarization), is shown in figure [5l Regions with 
plus (minus) signs depict up (down) spin phases for detector Note that the two 
regions are mutually exclusive so that while pure spin down is being detected in Di 
one can have either spin up or spin down in D2 depending on the range of a,p. The 
white regions correspond to no output at Comparing with figure Hb we see that 
no-output region are not identical for both detectors, these being larger for Di, i.e. one 
can have zero output at Di while having non-zero output at D2. As discussed before, the 
outputs depicted in figure |5] are also modulated by the magnitude of the corresponding 
component at the input, so the probability of the output exhibits noise coming from the 
random input spin orientation. 

5. Non diagonal mirror and beam splitter reflections 

Including the non diagonal matrix character of reflections at mirrors and beam splitters 
shifts the operation parameters of the spin filter but yields essentially the same 
qualitative results. The conditions must now be derived numerically. We start from 
equation m with the transmission and reflection matrices in equation |5l For the particular 
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Figure 6. The zeroes of the first (dashed line or red online) and second (solid line, 
blue online) eigenvalues oiUD2- When the first eigenvalue vanishes (and the second 
is non-zero), for specific combinations of a, /3 and ips the interferometer produces a 
perfectly polarized output in the |— ) state. Only a particular discrete set of solutions 
for is depicted. 

choice of 7r/4 incidence on the mirrors (see Appendix), the particularly simple non- 
tilting scenario described above is not possible. The extra parameter given by the angle 
of incidence on the mirrors/beam splitters lends itself to making this regime accessible, 
but we will not pursue it here. The more general scenario of a tilted axis yields a whole 
range of possible filtering solutions. 

Diagonalizing Wi52 in equation H] we find two eigenvalues. Setting the first eigenvalue 
to zero implies that in this rotated space the spinor is fully polarized (one of the entries 
of the output spinor is zero) as described in equation [HI Setting this eigenvalue to zero 
means setting its real and imaginary parts to zero. Such zeroes are depicted in figure 
|6]by the dashed lines (red online) for different values of the magnetic field and specific 
combinations A(a,/3), defined in equation |6l and 9 = tan~^(/3/a). In order for filtering 
to be performed such zeroes must be accompanied by non-zero values of the second 
eigenvalue in the same detector. The zeroes of the second eigenvalue are depicted in 
figure |6] by the solid lines (blue online) which are non-overlapping with the dashed lines 
for the first eigenvalue. Thus the figure shows alternative filtering conditions for either 
spin up or spin down in the tilted basis. 

The circular empty region in the middle of the plot correspond to non-polarized 
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Figure 7. The dashed curves represent zeroes of the first eigenvalue for tps = Stt/IOO 
upon a contourplot for the modulus of the second eigenvalue. The lighter shades 
represent higher values of the output polarization. One can extract the SO strengths 
from the plot by solving a simple system of equations for each value read off on the 
dashed curves. 

output in the tilted axis. Such a region contains some pointhke solutions that are of 
less interest experimentally since they would be difficult to tune. We recall that the 
previous discussion in section 14.21 is equally valid in this case, all incoming electrons 
at the input are polarized at the output no matter their energy as long as particular 
parameters ranges in the a, /3,^b space are met. 

In order to see the magnitude of the spin polarization for a particular value of 
the external magnetic field we draw a contour map of the magnitude of the second 
eigenvalue while the first one is zero. The background value at the dashed curves in 
figure [71 show the intensity of the pure down spin polarization at detector D2 when at 
ipB = Stt/IOO. The highest values of output achieved corresponds to the lighter shades 
on the contourmap. 

6. Summary 

We have proposed a perfect spin ffitering device based on a Mach-Zehnder type spin 
interferometer. The regimes of operation are subject to no limitations on the spin-orbit 
strengths and interferometer dimensions as in previous work. The treatment can be 
easily extended to unequal arm lengths and angles of incidence on the mirrors/beam 
splitters, that are likely to occur in the actual implementation of the interferometer. 
Such a generalization would provide additional parameters to manipulate ffitering 
conditions. In the simpler analysis above involving scalar mirrors, we find both a non- 
tilted and tilted axis spin filtering solutions referred to the axis of quantization in which 
one writes the input states and for arbitrary incoming energies. The non-tilted case 
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is not found in the scenario where the SU{2) gauge field is approximated by a U{1) 
hke gauge, and is pecuhar to the full non-Abelian treatment. This solution has the 
advantage of simplicity. On the other hand, the tilted axis solutions are shown to be 
well approximated by the Abelianized forms of reference |7j valid for certain reasonable 
conditions of SO strengths in relation to the interferometer arm lengths. When realistic 
mirrors/beam splitters are introduced, the mixing of the spinor components leads only to 
non-tilted solutions when 7r/4 refiections are contemplated. In this situation we run out 
of adjustable parameters to tune a non-tilted solution, that should be recovered when 
other incidence angles are considered. The qualitative scenarios for the operation of the 
diagonal and non-diagonal mirrors are the same and only the parameter combinations 
for filtering change. 

Perfect filtering means that all spins in one of the detectors are polarized always 
in the same axis and orientation. This has the drawback that the current is not steady 
since the probability of producing a completely polarized electron varies with the initial 
projection, of the input spinor, onto the chosen quantization axis. This projection is 
random as electrons are injected from the Fermi seafSOj. A density matrix approach 
should be implemented so that one can also assess finite temperature effects on the 
filter operation. It should be also noted that the interference setup does not produce a 
pure spin current, since polarization is accompanied by a charge current. 

An interesting insight, exploiting the analogy with the Aharonov-Bohm effect in 
the Abelian case, comes from observing the role of A in the non-Abelian case. A and the 
voltage V essentially play the same role as the pair and magnetic flux. Indeed, 

for a purely Pauli type SO interaction, as A = {mL/h)a and a = heE/{m^c^), then A 
can be rewritten as 27ri?L/(27rmc^/e) = 2t[Ve/Vq, where Ve = EL, the voltage along 
the arm of length L in an electric fleld of strength E. Vq is a quantum of volt age p^. 
Although Vq is very large for this calculation, the material Rashba coefficient would 
lower it to the order of 1 eV/e. 
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Appendix 

Here we derive the general conditions for refiection at a beam splitter on a mirror in 
the presence of both Rashba and Dresselhaus interactions. Starting from Hamiltonian 
in equation [T] we can solve exactly for the eigenvalues and eigenfunctions. Ignoring the 
Zeeman term we have 




(16) 
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Figure 8. Detector D2 output while Di filters out the up spin component (spin down 
polarization) . 



with eigenfunctions given by 



|k ±),, 




;i7) 



-fF'^kx, ky 

where k = {k^, ky), ± stand for the two eigenvalues and the subindex i stands for incident 
wave. The convention we take, according to the figure, is that k^ and ky are positive 
components for the incident electron. Referred to those components, one can obtain the 
reflected basis components by changing k^ — ?■ —k^ and ky — )■ ky as the momentum in 
the y direction is conserved. To obtain the projections in terms of the reflected basis we 
write 

|k ±)i = a±|k +)^ + 6±|k -)^, (18) 

where the subindex on the right indicates the reflected complete basis set. One can 
then compute the superposition coefficients a± and b± by performing the appropriate 
overlaps between incoming and outgoing wavefunctions 



a± 
b± 



^{k + \k±), = 1/2 [1 ± F*{-kx, ky)F{kx, ky)] , 

^{k - \k±), = 1/2 [1 T F*{-kx, ky)F{kx, ky)] . (19) 

Each of the outgoing amplitudes gets multiplied by the scalar reflection coefficient r in 
the case of the beam splitter and r = 1 for perfect mirrors. The previous coefficients 
govern the QPCi, the upper reflection of QPC2 and Mi in figure [H while exchanges 
of kx — !■ —kx would generate the corresponding matrix for the M2 and the bottom 
reflection of QPC2. 

The wavector components can be expressed as k = (A; sin 7, /c cos 7) for a generic 
incident angle as seen in the flgure. For the case of 7 = 7r/4, the reflection matrices 
are particularly simple and one obtains equation [5l where the transmission matrix is 
trivially diagonal since the electron beam does not change direction. 
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A coordinate independent way to state the general result is by identifying 

F{kx,ky) = expi(j)i and F{—k^,ky) = expi(j)r then one can write the full 
reflection/transmission matrix as 

/ rcos[(0^ - (/)i)/2] ir sm[{(f)r - /2] t \ 



ir sin[(0r — 0j)/2] r cos[(0r. — t 

t r cos[(0r — 0i)/2] — ir sin[(0j, — 0j)/2] 

\ Q t -irsin[(0^ - 0i)/2] r cos[(0r - 0i)/2] I 
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